Entanglement and control operations in Ising interactions of bipartite qubits 
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Entanglement generated by Ising model has been studied for several authors in order to understand 
the relation between it and magnetic properties of materials, principally using one or two dimensional 
models for two or more particles. In this work, Ising model evolution is solved in three dimensions for 
two parts including an inhomogeneous magnetic field, giving an detailed study of the entanglement 
properties derived from interaction. Some relations between entanglement and energy or spin are 
developed specifically to stablish candidates for entanglement witness. Finally, some basic quantum 
control operations are prescripted for this model as to preserve the properties of the system as for 
transfer information between the two parts. These last schemes of control are useful when initial 
state is completely known, to gain domain on the system. It prevents to introduce more sophisticated 
control schemes which normally are necessary when initial state is at least partially unknown. 

PACS numbers: 03.67.Bg; 03.65.Ud; 03.67.-a 



I. INTRODUCTION 

Entanglement is used in quantum computation as a 
central aspect for improve information processing in or- 
der to exploding the interesting properties of quantum 
mechanics [l|, [1]. For this reason, entanglement is sub- 
jected to deep research to understand completely its com- 
plexity, properties and potential usefulness 3, 4, 5]. In 
the last sense, its control is one of the most important 
aspects of it p, [zl , nevertheless its study will not have a 
complete map of road until his quantification and behav- 
ior could be understood. 

Nielsen Q was the first to report studies of entangle- 
ment in magnetic systems based on a two spin systems 
using the Ising model with an external magnetic field. Af- 
ter of this, different studies have extended this research 
for more complex systems and depending on external pa- 
rameters (as temperature and strength of external field) 
^, 10], and considering different models of Ising interac- 
tion (XX, XY, XYZ depending on focus given by each 
author in order to reproduce calculations related with 
one or two dimensional lattices) [ll|, [H, [H, [H, . 

In this paper we just study bipartite systems in three 
dimensions, but more in the aim of learn about entangle- 
ment and control related with Ising interaction and infor- 
mation transfer processes, than to study lattice proper- 
ties of entanglement. Extensions to multipartite systems 
will be studied in future works. 



II. ISING INTERACTION AND VARIANTS 

Ising model is motivated mainly by far-field strength 
of a magnetic dipole interaction between two particles in 



which the energy of bounding is given by: 



E 



(nil • m2 — Siiii • fm2 • f ) 



(1) 



where r is the vector distance between particles, r its 
unitary vector associated and xn.i is the magnetic mo- 
mentum of particle i. In addition, a model which relates 
the magnetic momentum with the spin is: 



m, 



(2) 



with gi as a tensor. Combining two last expressions we 
obtain: 



E 
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Depending on g^ij values the Ising-like interaction be- 
come in different interaction models which have been 
proved by several authors in studies of entanglement 
fill . m, [H, [H, . Ising model is precisely which where 
Jji is proportional to the identity (see appendix A): 



E = — JSi • S2 



(4) 
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This kind of interactions were first used in statisti- 
cal physics to describe the magnetic behavior of lattices 
in different ways precisely by Ising ^Td, ^V?\ and after 
by Heisenberg [18j | in quantum mechanics introducing 
Pauli matrices proportional to f^. Works listed before 
(Til IT^ . [isl . [13, [iBl and other more recent showing trans- 
ference and control of entanglement in bipartite qubits 
[l9l | and lattices [13, normally focus on describe mag- 
netic properties of materials related with entanglement. 
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By other hand, Cai [24I has considered a more general 
model with Jji diagonal in order to study the relation 

between entanglement and local information. Closer this Ei = — J — 

last, but related with the previous papers, this work is = —J-\-B^ 

focused on the full Ising model in three dimensions, in E — J — R 

order to explore some properties of entanglement which 

arise between particles inclusively by adding an interac- 1^4 — j -t -n, \^i>) 

tion with an inhomogeneous external magnetic field in 

the z direction: ^^^^ eigenvectors: 



H 



-J<Jl ■ (72 + Si (712 + B2a2:^ 



(5) 



and some operations that one could apply to control en- 
tanglement. 

Using explicit form of Pauli matrices we obtain in ma- 
trix form for ([S]): 



H = 



B+-J 
B_ + J 
-2J 





(6) 



where B+ = Bi + B2, B_= Bi- B2, R = Bl + 4J2. 

One more suitable selection of parameters is (to make 
finite some parameters and to reduce some expressions): 



b+ = B+/R,b^ =B^/Re [-1,1], 
j = J/Re [0,l/2],t' ^ Rt 



(7) 



We will not drop the prime in the time because the 
letters should be advertise to the reader about the actual 
selection of parameters. We use it in all cases which it is 
situable, except for some final expressions. 

Diagonalizing the Hamiltonian ([5]), we can obtain the 
operator of evolution, which become in Dirac notation: 

U{t') = e-'*'(''+-^') IO1O2) (O1O2I 

-He-'*'^ (cost' - ib- sint') IO1I2) (O1I2I 
-Hi2je"**'j sint' IO1I2) (I1O2I 
+i2je-"'^ sint'|li02) (O1I2I 
+e-'^'^{coiit' + ib^ sini') II1O2) (I1O2I 
+e'*'(''++-'")|lil2)(lil2| 

(8) 



III. EVOLUTION AND PROPERTIES IN THE 
ISING MODEL FOR 2-PARTS 

A. Generalities 

It is well known that Ising-like interaction generate en- 
tanglement (Tlj . Our Hamiltonian has the eigenvalues 
(Fig. 1): 
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(10) 



because of that, these states are invariant under interac- 
tion Ml. 



B. Energy, spin and entanglement 

Spectcd values of energy and spin are normally related 
with entanglement. Note first that just for some selec- 
tion of parameters (when ^ j ^ 1/2), the two 
last states of ([T0| will correspond to maximally entan- 
gled states: IM3) = |/3oi) and \u4) = |/3io), it means in 
an homogeneous field, this Bell states in z direction ba- 
sis are invariant. Present section present some remarks 
about relation between energy and spin with entangle- 
ment. We will begin since the most simple until most 
general scenarios. 

Firstly, in absence of field one can show by direct cal- 
culation that for 1^) separable, then (see appendix B): 



E 



\H\^)^ [-J,J] 



(11) 



it means that Hamiltonian could be used to define an en- 
tanglement witness operator because of (/3io |i?|/3io) — 
3J. Effectively, in this case IO1O2), II1I2) (or alterna- 
tively, |/3oo)j |/3oi), form an invariant subespace of 
TC^"^ (all state obtained by combine this states is inva- 
riant) isolated from |/3io). In this case, energy is degene- 
rated for any state of three first and it is £'1,2,3 = —J^ 
instead for the fourth state it's £4 = 3 J. In absence of 
field, non invariant evolution it's only possible by com- 
bining both subspaces. 

After, by adding an homogeneous field, two of the three 
first degenerated cigenstates become unfolded with ener- 
gies El = — J — B+, E2 = —J + 5+ (depend on mag- 
netic field, B — i?+/2) with separable eigenstates IO1O2), 
II1I2) respectively. Instead £3 = —J, £4 = 3J remain 
unchanged and independent of magnetic field with the 
same entangled eigenstates |/3oi) and |/3io), which are now 
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the only invariant Bell states in z direction basis. This 
apparent difference between influence or non influence 
on field appear as reason to get or not to get invariant 
separable eigenstates under evolution operator ([5]). 

In the general case when magnetic field is inhomoge- 
neous, energy of the two separable eigenstates unfolded, 
remain now unchanged. Nevertheless, the previous maxi- 
mally entangled eigenstates become just partially entan- 
gled with energies E-^ = J -R< - J, Ei = J + R>'iJ. 
Otherwise, one can prove than the two last eigenvectors, 
kta) and |u4), are not separable at least that — s- oo 
[27| . By direct calculation one can show that for j?/;) 
separable (see appendix B): 



E ^ {iP\H\il;) (^[-M,M] (12) 

where M = Max{ J + |-B-|, In this case, energy is 

no more in all cases an entanglement witness because for 
(/3i^o| H |/3i,o) = 3J B_ = EiQ. With this, one of M, 
-E4 and -Eio could be greater values than other, depen- 
ding on values of i?_ and i?+ . In some sense, complexity 
and difference in intensity of the magnetic field appear 
to destroy this relation energy-entanglement. 
By comparison, the eigenvalues of 3i ■ ai are: 



|0l02)(0l02|-f 

-(l-f8j6_ sin^t') |0il2> (O1I2I + 

2(1 - 26^ sin^ t' sin 2i') IO1I2) (I1O2I + 

2(1 - 2\?L sin^t' - sin2t') II1O2) (O1I2I + 

-(1 + 8j6_ sin^ i!) II1O2) (I1O2I + 

|lil2> (I1I2I 

(17) 

then for energy eigenstates ([T(I| -.ax-a^ become still time 
independent: 

(Ui |(Ti • (T2I Wl) = 1 
(m2 lo'l • i?2| W2) = 1 

(w3 |cti • CT2IW3) = 4j-le(-l,l] 
(w4 |(Ti • 1T2I U4) = -4j - 1 e [-3, -1) 

(18) 

noting that the last eigenstate gives values corresponding 
to non separable states, so this quantity is a real entan- 
glement witness. 



Al,2,3 = 1,A4 = -3 
with the eigenvectors: 



W) = IO1O2) 

^2) = I1112) 
1 



(13) 



\v^) = -^(|0il2) + II1O2)) = |/?oi) 
W) = ^(|0il2)-|li02)) = |/3io) (14) 



An important issue is that if some state jV') is separa- 
ble, then: 



(cti -02) = (V^l^i •^2|V') e [-1,1] 



(15) 



so (7i • (72 can be used to construct an entanglement wit- 
ness operator in this case because (because of his propor- 
tional dependence with energy seeing before for abscence 
of field case) : 



(cti • $2) = (/3oi| CTi • CT2 |/3oi) = -3 



(16) 



The evolution ^ of {ai ■ 02) become dependent just 
of B_: 



C. Entanglement and separability 

Using evolution operator we can verify that neverthe- 
less IO1O2), II1I2) are invariant, the states: 

C/(i')|0il2) = e-^^'*'((cost'-i6_sini')|0il2) + 
2ijsint' II1O2)) 

[/(t')|li02) = e-'^*'((cost' + i6_sini')|li02) + 
2ijsint' IO1I2)) 

(19) 

have an interesting behavior. By calculate Schmidt coe- 
fficients, we note that these states are maximally entan- 
gled when: 



cos^ i! + }?_ sin^ il = Af sin^ t' ^ B't = -4 cos 2Rt 



(20) 



this last, for non normalized time. It means: 



ta 
tb 



1 / Bl 

— arccos ^ 

2R \ 4J2 



nT 
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R 2R 



■ arccos 



'4J2 ] 



-|- TzT, n e 



(21) 
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\2§\. This happens 



the period of the process is T 

only if < 1. Resulting states after of these times 
of evolution are respectively (dropping some unitary fac- 
tors): 



C/(ia,fc)|0ll2) = 



i , iarctan(-i^-g sgn(tan _Rt„_b)) |n i \ , h n \^ 

-^(e - IU1I2) + II1U2)) 

U{ta,b)\ll02) = 



r , — zarctan(-^^ — ^ sgn(tan 



72 



(e 



I1O2) + IO1I2)) 
(22) 



where sgn(a;) = if a; 7^ 0. They reach two times 
maximal entanglement for each cycle, except if = 
2 J. In this case, the last both states become |/3oi), |/3io) 
(just one of them depending on sign of B_) at different 
times periodically. Inversely, this is the condition for the 
last Bell states become separable. 

Figure 2 shows the Von Neumann entropy for diffe- 
rent initial states under this evolution (graphs are equal 
for IO1I2) and for II1O2), nevertheless different states are 
reached) . 

We don't show here evolution for initial maximally 
entangled states generated combining IO1I2) and II1O2), 
but the results are similar. When \B^\ — 2 J this state 
reaches with the same period a separable state, other- 
wise just reaches partial entanglement state. In the same 
way, they are invariant when £?_ — > or periodical when 
B_ 00. 

Its important remark that the periodic behavior shown 
by last initial states it's only a partial view of phe- 
nomenon. Because of different energy eigenvalues arosen 
with the inhomogeneous magnetic field don't give ratio- 
nal quotients necessarily, this behavior for states initially 
entangled or separable have a non-periodic evolution in 
general. By example, taking as initial state ( 9 parameter 
goes from to ^ , taking maximal entanglement in both, 
but having initial partial entanglement in intermediate 
values): 



i/j) = sin9 |/3oi) - cose |/?io) 



(23) 



and let it to have Ising interaction, we find that the 
Schmidt coefhcients become after time t': 



Ai,2 = -(I±(16j2(l-4j2)sin4t'sin4 0-f 

sin^ 20(sin^ 2jt' + sin^ t' cosAjt' - j sin4ji'))^) 

(24) 



Here, |&-| = y/l—Tp has been used. Figure 3 show 
entanglement evolution of this more interesting state for 



different values of 6 and j, which exhibits the properties 
of this interaction. Appendix C shows some details about 
properties of entanglement evolution and conditions of 
periodicity of the separability for general cases. 

By comparison with results of the before subsection, 
one can calculate using (fTT]) for IO1I2), II1O2): 



(0ll2|t?l •c?2(i')|0il2) 

(ll02|c?l •C?2(O|ll02) 



-{l + 8jb-sin^t) e [-3,1] 
-(l-8.?6_sin2t) e [-3,1] 
(25) 



because 8j5_ e [—2,2]. {ai ■ (J2{t)) = —3 when 8jfe- 



in t = 



_ 2rt+l 



_ 2k+1 



2j^,n G "L for the first state and 8j6_ = —2 in 
2^ , n e Z for the second state. These are precisely 



the times when these states become |/3io)- In virtue of 
our previous result about entanglement witness, it tell 
us that this initially separable states reach an entangled 
stage by the action of the magnetic field . Precisely this 
is the reason because of one invariant state in the absence 
of magnetic field, now with inhomogeneous magneticfield 
we can achieve an separable form. 



IV. ELEMENTAL PROCEDURES OF CONTROL 
A. Evolution loops 

Evolution loops were introduced by Mielnik [11] and 
applied and extended in other directions of control by 
several authors [l^, [H, as simple operations to pur- 
suit the specific behavior of quantum system. Note that 
here we don't introduce any stochastic element as con- 
sidered in @, 0] ■ 

In some different sense we can reproduce this kind of 
effects by impose that the evolution operator after of 
some time T of application become: 



U{T) 



I 



(26) 



where / is the identity operator. Taking the evolution 
operator ([5]), we note that the first condition that needs 
fulfill is RT = mr. In addition, analyzing the diagonal 
terms to fit to the form (|26p . we need impose the condi- 
tions: 



B. 



2{s - m)J 
m + s — 
2J ^ {2n — m — s){m + s) 



B_ = ±- 

m + s — n 

{m + s — n)'K 



T 



2J 



(27) 



with n, TO, s G Z and 0<?i<to-|-s< 2n. We obtain: 



U{T) = (-l)"e 



n-i JT 7 



(28) 
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It's clear that ii p G Z+, then all cases with: m,s,n ^ 
pm,ps,pn are the same process but with T — > pT. Other 
cases are not physically equivalent. The fastest process, 
with m + s — n — 1, appear as requiring stronger fields, 
nevertheless still is available without magnetic field by 
choosing n,m,s = 1 : B+ = B- = 0,T = jj. Figure 4 
shows the entanglement evolution for the family of states: 



l</'±) =\^|0il2)± Vl^|li02>,pe [0,1] (29) 
under these evolution loops. 

B. Transference of information 

Another control operation which could be induced is 
the transference or exchange of information between par- 
ticles. Suppose that two particles in separate states begin 
to interact trough Ising interaction (O. It's possible that 
after some time T this particles exchange their states? 
In order to get this effect, the evolution operator should 
be: 



U{T) = e-'^^h^^ 



(30) 



where is a phase and /i^2 is the unitary exchange 
operator between particles 1 and 2: 



/l^2 = |0l02)(0i02| + |0il2)(ll02| 
II1O2) (O1I2I + II1I2XI1I2I 



(31) 



This effect could, by example, induce the information 
transference: 



for the fastest process with n = 0. This last it's only 
the fact that repeating two times this process one obtain 
of course an evolution loop, so this kind of operation is 
a subclass of Evolution loops but with one half of its 
period. 

Figure 5 exhibit time evolution of entanglement for di- 
fferent initial separable states in one period of this effect 
depending on p and t. It's remarkable that when states 
are less similar (more information to exchange), interme- 
diate entanglement have stronger variations. Otherwise 
in the case of initial separable state the process require 
intermediate entangled states for the exchanging. All 
it means that while more information it's needed to ex- 
change, systems require more intermediate entanglement 
increase. 



V. CONCLUSIONS 

Study of the interaction between single pair of qubits 
is rich in complexity, specially when magnetic field is in- 
troduced. Erratical behavior is one of the principal as- 
pects of this complexity because generalization to larger 
Ising strings is not easy, specifically that related with en- 
tanglement. Properties of the system exploited trough 
quantum control could let to emerge some useful aspects 
related not only with his driven but with transference of 
information, at least in the most direct approach (alter- 
native or maybe related with presented by p^ ) because 
still the non-linear aspects of these strings could keep 
much more secrets and surprises. Study of control to 
stabilize this kind of systems and to drive it in order to 
transfer information in larger chains or rings could be 
considered trough of operations presented here. 



U{T){a\Oi)+p\U))(^\i;), 



e*^|V)i<5§(a|02)+/3|l2)) 

(32) 



Fitting the Ising evolution operator to this last one, we 
obtain the conditions: 



B+ = 

B_ = 

T = 



2B, 
0, 

{2n- 



8Jm 
2n + l' 



l)7T 



4J 



(33) 



with n e Z+, 7T!, e Z so the evolution operator becomes: 



uiT) = ti-iy 



-iJT 



h^2 



(34) 



It's remarkable that this phenomenon happens just 
with homogeneous field. In addition this effect could hap- 
pen without magnetic field (choosing to = 0), inclusively 
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APPENDIX 

Appendix A: Symmetric Ising-Heisenberg model 
reliability 

A correct selection of parameters gtij shows that @ 
it's always possible. By example selecting a symmetrical 
relative position of particles: f = ^(1,1,1) and: 
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= -522(.9ll + 531) - .912(521 + 531) - (ffll + 52l).932 

-.923(.9ii + .931) - .913(521 + 531) - (511 + 521)533 
= -523(512 + .932) - 5i3(.922 + 532) - (512 + 522)533 

512 = -C,5l3 = C,g23 = -C 
521 = C,531 = -C,532 = C 

(35) 

with c constant. If c G M then J > 0, and if c G I then 
J < 0. 



Appendix B: Range for spected value of an 
observable for bipartite separable states 



/3(i 



(/3(cost' - ib^ sint') + 2ijjsmt') 



lit') = 

g-ut'(-^2ij sint' + 7(cost' + ib- sint') + ■j2ijsmt') 
5{t') = 5e*(''++^')*' (41) 



the eigenvalues of: 



a b 

c 



wiU be the Schmidt coefficients: 



(42) 



The problem of evaluate the range of an observable 
for a bipartite separable state reduces to consider the 
function: 



A,,, = -(l±Vl-4|A(t')|4) 



with 2£ 



(43) 



f{x,y) = 2(x^l-y^+y^Jl-x^f + 

I3i{2x' ~ I) + I32{2y^ - I) (36) 

with x,y [0, 1] and solving an optimization problem 
using calculus, which conduces to solutions for the ex- 
trema: 



{-/3+,2-/3_,2 + /3_,/3+} (37) 
- ^i±^ and (3- = ^i^, when a; = 0, 1 



where (3+ — ^ 
and/or y — 0,1. 



|A(t')p = {ac-\b\') = \a{mt')-m7{tT (44) 

\ip{t')) will be separable at time t' iff A(i') = 0. After 
some calculations we obtain: 



A(t') = aSe^'^'' - (3je-^'^'' - 

{13^ - 7^)2j6_ sin^ t' - ^p-ff sin^ t') (45) 
Using A = A(0), /3 = 7re"^ one can write (if 7 7^ 0) 



Appendix C: Periodicity of entanglement and 
separability in Ising model for two qubits 

Taking a general bipartite state: 



1(^(0)) = a IO1O2) + (3 IO1I2) + 7 II1O2) + <5 II1I2) (38) 



with |ap 
in matrix form as: 



P + |7p + |(5p = l. Representing this state 



A[Q) 



a (3 
7 5 



(39) 



applying the Ising evolution we obtain for it: 



A{t) = 



a{t) Pit) 
lit) S{t) 



(40) 



where: 



F{r,cl),j,t') 



i{<t>-<t>A 



-(A(t')~ Ae^'^*') = 



7" 

-2r^ije"t'-^^^*' sint'(cost' - ib-{j) sint') + 
2r{i sm2jt' + ie^^'^*' f sin^ t) - 
2ije-"^-^'^^' sint'(cost' + ib_{j) sint') 



(46) 



where b^{j) denotes the dependence of &_ from j. 

Here (f)^, (f>A' are the phases of A, A(t') respectively. 
This shows the non periodic behavior of separability and 
entanglement at least that j e Q. When the right side 
of this equation become zero for t' > and initial state 
is separable, then we have separability at time f again. 
But from the last formula, the same is true for any other 
value of |A|, because if F{r,(j),j,t') vanishes, it implies 
that the state (in general, different of initial state) at 
time t' have the same value: |A(t')| = |A|. In addition 
their relative phase should be equal to 2jt' until some 
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multiple of 27r. Figure 6 show behavior of some cases of 

F{r,cl),j,t'). 

Some properties of F(r, (j), j,t') are remarkable. If 
F{r,(j),j,t') = Fr{r,(j),j,t') + iFi{r,(j),j,t') are the real 
and imaginary parts of F, then for some r and (j) fix, we 
have that if F{j = x,t' = t') = it means that until 
first order: 



maximally entangled state, precisely |/3oi) or |/3io)- Note 
too that while magnetic field is more inhomogeneous the 
process is faster but far away of limit case \B-/J\ = 2, 
the maximum of entropy reached goes to zero. 



Figure 3 



J^r(X + dj,T'+dO = 

Fi{x + dj,T' + dt') 



dFr 

m 

dj 



OF 

m 

dt' 



dt' 



dt' 



Entanglement evolution for different values of 6 and 
J = ^ € [0, i]. Periodic behavior is only possible for 
rational j: a)jQ, b)^, c)\ and d)^. Increasing j we in- 
crease oscillations between values of entropy of entangle- 
ment, while for 9 = ^ the maximum variation is reached 
(47) in all cases. Evolution for j irrational is shown in e), with 



If Frix + d.7, t' + dt') = F,{x + dj, t' + dt') = still. 



the condition of roots preservation, 4^=0, gives: 



dFr 




dFi 




dFr 




dFi 






dt' 




w 




dj 



= 



X,t' 



(48) 



In our case we obtain the result for last equation (ta- 
king t' = mn, X = ^) the cases where we know that F 
is zero for rational solutions): 



J = Ri I , showing a non-periodic behavior. 



Figure 4 

Time evolution of entanglement under evolution loops 
in one time period for different values of p. a) The case 
m,n,s — 1 for |i^±); the case n = 3,m = 0, s = 2 ior b) 
\ip+) and c) \(p-)- Note in this last case the equivalence 
of inversion on the p axis. In spite of that there are more 
than on cycle of periodicity of the graphs, the original 
states remain unchanged just when t = T. 



- 47rnr(r^ - 1) sin ^ = (49) 

it means that for cases r = O,lor0 = O,7r the rational 
solutions t' = rmT,j = ^,n,m € Z are preserved (not 
necessarily with the same values) for j € Q' (see d-f in 
figure 6). 

FIGURE CAPTIONS 
Figure 1 

Eigenvalues of energy depending on magnetic field 
{B-/J and parameters). Since Ei/J (darkest) 

until E4/J (lightest). Note how E4/J is separated, sug- 
gesting different behavior which can be verified trough 
some properties. 

Figure 2 



Figure 5 

Time evolution of entanglement under exchange of in- 
formation depending on p-parameter and time t, using 
■m,n — 0. For the initial state: a) \(pa) = {y/p\Oi) + 

VT^IW) ® I02), b) = (v^ioi) + yr^lii)) 
(Vp|02)-yr^|i2», c) 1^,) = (^/^^|Ol)+vp|ll))® 
(Vp|02)-yr^li2)), d) = (VT^|Oi)+vp|ii»0 

{y/p IO2) + \/l — p |l2))- Graphs suggest that when more 
information should be transferred between parts, then 
more entangled become in the intermediate stage. Fi- 
nally, just for comparison, the evolution for the entan- 
gled states: e) |ye) = v ^|0i02) + VT^lhU), f) 
\(pf) = ^|0il2) + \/l — p II1O2). In the first case en- 
tropy of entanglement is unchanged, but in the second 
slightly variations appears when partial entanglement is 
weaker doubted to the same exchange of information- 
entanglement phenomena. 



Entropy of entanglement behavior for the evolution of 
IO1I2) and II1O2) for different values of \B-/J\ in which 
can be reached maximally entangled states. Using J = 1, 
periods shown are different: a) \B-/J\ = T = tt, 
b) \B./J\ = 1 ^ T = f^, c) \B./J\ = 2 ^ T = 

d) \B_/J\ = 4 T = ^. For higher values, as 

in d), just partially entangled states are reached. Note 
in c), the limit case, that just two times is reached a 



Figure 6 

Behavior of \F{r,4',j,t')\ for a set of values of r (de- 
creasing with darkness): a-c) shown for rational j = 1/4 
exhibiting 27r periodicity of function and roots; note that 
(j) is related with symmetry between roots, d-f) for irra- 
tional j, where periodicity is lost, it is remarkable that 
some roots appear accumulated near from t' = Anp, p € Z 
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(specially for the cases r = 0, 1 or = 0, tt in accor- 
dance with result ([H])). For rational j it corresponds 
with m — 4, which implies j — j^^- Effectively, taking 

n = 3 we note that j = | 0.375 w 0.377 = -j=. Any- 



way, repeatability of values for entanglement is normally 
present in Ising model, but they don't follow a periodic 
pattern. 
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